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Abstract: 

Decay rates of quarkonia are studied within the framework of NRQCD formalism. The basic parameters 
of the formalism have been obtained from different potential schemes studied for the spectra of quarkonia. 
We estimate the heavy quarkonia mass spectra, radiative and leptonic widths and compare them with other 
contemporary theoretical approaches and experimental results. 



1 Introduction 



The spectroscopy and decay rates of quarkonia are quite important to study due to huge amount of high 
precession data acquired using number of experimental facilities at BES at the Beijing Electron Positron 
Collider (BEPC), E835 at Fermilab and CLEO at the Cornell Electron Storage Ring (CESR), the B-meson 
factories, BaBar at PEP-II, Belle at KEKB, the CDF and DO experiments at Fermilab, the Selex experiment 
at Fermilab, ZEUS and HI at DESY, PHENIX and STAR at RHIC, and NA60 at CERN pQ. New states and 
production mechanisms, new decays and transitions have been identified and even larger data samples are 
expected to come from the BES-III upgraded experiment, while the B factories and the Fermilab Tevatron 
will continue to supply valuable data for few years pQ. New facilities like the LHC experiments at CERN, 
Panda at GSI etc. will offer greater challenges and opportunities in this field [HO [3]. 

The mesonic states are not only identified with their masses but also with their various decay rates. So, 
one of the tests for the success of any theoretical model for mesons is the correct prediction of their decay 
rates. Many phenomenological models predict the masses correctly but overestimate the decay rates [4j 
G3 [51 [7J [HI M HOI EI]- For better estimates of their decay rates with reference to the experimental values, 
various corrections due to radiative processes, higher order QCD contributions etc have been suggested [12] , 
In this context, the NRQCD formalism is found to provide systematic treatment of the perturbative and 
non-perturbative components of QCD at the hadronic scale [321 [HI US] . For the present study, we employ 
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phcnomcnological potential schemes for the bound states of heavy quarkonia and the resulting parameters 
and wave functions to study the decay properties. For example, we include here the extended relativistic 
harmonic confinement scheme |161 117] employed for computation of quarkonia masses and decay rates along 
with other potential approaches such as coulomb plus power potential (CPP„) with different choices of the 
power index v etc [9] to study the meson properties. 

2 Potential model approaches for quarkonia spectroscopy 

In this section, we review important features of the models such as ERHM and coulomb plus power poten- 
tial (CPP„) employed in the spectroscopic study of heavy flavour hadrons. We extract the spectroscopic 
parameters of the models such as the model parameters that reproduces the quarkonium spectra, the wave 
function at the zero separation etc. 

2.1 Extended harmonic confinement model for hadron spectroscopy 

Choice of scalar plus vector potential for the quark confinement has been successful in the predictions of 
the low lying hadronic properties in the relativistic schemes for the quark confinement [HI [19l [20 . In 
relativistic harmonic confinement model (RHM), coloured quarks in a hadron are confined through the 
action of a Lorentz scalar plus a vector harmonic oscillator potential. It also leads to an alternate scheme 
for the understanding of nucleon-nucleon interactions from a more fundamental level [21) . The RHM has 
been extended to accommodate multiquark states from lighter to heavier flavour sectors with unequal quark 
masses [TBI IT?] . 

The mass of a hadron having p number of quarks in this extended RHM (ERHM) is expressed as [THl [T7] , 

p p 
M N (qiq 2 ) = },e N (qi,p) con f + 2J ^N(qiqj)coui 

i—l i<j—l 
P 

+ e ivfe>?i)sD (1) 

i<j=l 

where the first sum is the total confined energies of the constituting quarks of the hadron, the second sum 
corresponds to the residual colour coulomb interaction energy between the confined quarks and the third 
sum is due to spin dependent terms. 
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The intrinsic energies of the quarks using harmonic potential is given by 



\ 



3M 

(2N + 3)n N (q) + M\ - / q n ( g ) (2) 



E 

i=l 



M g (i) 



The coulombic part of the energy is computed using the residual coulomb potential given by 22J, 

V cou i(q i qj) = (3) 

UJ n T 

here io n represents the colour dielectric "coefficient" [22 which is found to be state dependent [T7J, so as 
to get consistent coulombic contribution to the excited states of the hadrons. Such state dependence in 
the effective a s for the excited states of quarkonia has been reported by others 23J . It is a measure of the 
confinement strength through the non-perturbative contributions to the confinement scale at the respective 
threshold energy of the quark-antiquark excitations. 

The wave functions for quarkonia are constructed here by retaining the nature of single particle wave function 
but with a two particle size parameter fl^(qiqj) instead of &,N(q) [16] and the coulomb energy is computed 
perturbatively using the confinement basis with the two particle size parameter. The spin average (center of 
weight) masses of the cc and bb ground states are obtained by choosing the model parameters m c — 1.428 GeV, 
m b = 4.637 GeV, k = 0.19252 and the confinement parameter A = 0.0685 GeV 3 / 2 [17] . 

From the centre of weight masses, the pseudoscalar and vector mesonic masses are computed by incorpo- 
rating the residual two body chromomagnetic interaction through the spin-dependent term of the COGEP 
perturbatively as, 

4i{*qj)s.D. = (NJ\V SD \NJ) (4) 

where \NJ) is the given hadronic state. For mesons \NJ) becomes the \qiq~j) states. We consider the two 
body spin-hyperfine interaction and the spin-orbit interaction of the residual (effective) confined one gluon 
exchange potential (COGEP) given by [24], 

a s (fi)N?Nj X l .\ J 

4 [Ei + mi}[Ej +m j ] 

x[^6 3 ( nj ) - C A OOM r 2 D 1 {r lj )]{- 2 -a i ■ a,) (5) 



Va 



and 



v ls N i N ! A, -A. 



q ' q > 4 (Ei + MiXEj+Mj) 2r l3 

x iL ■ S (D' (r ij ) + 2D' 1 (r«))l (6) 
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Table 1: Masses of charmonia in GeV/c 2 



State 


ERHM 




CPP„ 




32 


M 


M 


[35] 


M 


[37J 






i/=0.5 


i/— 1.0 


v— 1.5 


i/=2.0 


PDG 












l L S 


2.985 


3.000 


2.950 


2.912 


2.882 


2.980 


2.981 


2.990 


2.979 


3.093 


2.982 


2 l S 


3.626 


3.352 


3.522 


3.636 


3.852 


3.638 


3.625 


3.627 


3.588 


3.096 


3.619 


3^0 


4.047 


3.541 


3.912 


4.212 


4.436 


_ 


4.032 


_ 


3.991 


_ 


4.053 


4 1 S 


4.424 


- 






- 


- 


4.364 


- 


- 


- 


- 


l^i 


3.096 


3.092 


3.112 


3.129 


3.144 


3.097 


3.089 


3.097 


3.096 


3.096 


3.097 


2 3 S 1 


3.690 


3.375 


3.583 


3.739 


3.852 


3.686 


3.666 


3.685 


3.686 


3.476 


3.686 


3 3 S 1 


4.082 


3.553 


3.950 


4.285 


4.547 


4.040 


4.060 


4.050 


4.088 


3.851 


4.102 


4 3 5i 


4.408 










4.415 


4.386 


3.443 




4.223 


4.447 


1 3 P 


3.431 


3.329 


3.398 


3.461 


3.479 


3.415 


3.425 


3.496 


3.424 


3.468 


3.415 


l 3 Pi 


3.464 


3.302 


3.424 


3.504 


3.556 


3.511 


3.505 


3.525 


3.510 


3.468 


3.511 


1 3 P 2 


3.530 


3.323 


3.477 


3.590 


3.673 


3.556 


3.556 


3.507 


3.556 


3.467 


3.556 


l 1 ^ 


3.497 


3.313 


3.450 


3.547 


3.615 


3.526 


3.524 




3.526 


3.467 


3.524 


2 3 P 


3.891 


3.494 


3.792 


4.009 


4.153 


3.800 


3.851 




3.854 


3.814 


3.864 


2 3 Pi 


3.899 


3.507 


3.832 


4.085 


4.269 


3.880 


3.923 




3.929 


3.815 


3.950 


2 3 P 2 


3.916 


3.531 


3.911 


4.239 


4.501 


3.940 


3.970 




3.972 


3.815 


3.992 


2 x Pi 


3.907 


3.519 


3.872 


4.162 


4.385 




3.941 




3.945 


3.815 


3.963 



Where D' (rij) and D[(rij) appeared in Eqns ([5]) and ([6]) are derivatives of the confined gluon propagators 
of CCM [53] and Ni/j are the RHM normalisation constants given by, 

Ni/j = V /2(P 4/J + M lM )/{iE l/3 + M lM ) (7) 

Where, C CCM corresponds to the confinement strength of the gluons and r y is the interquark distance. The 
computed masses for charmonia and bottomonia based on this ERHM model are given in comparison with 
experimental and other theoretical model results in Tables [T] and [5] respectively. The resultant spectroscopic 
masses and their radial wave functions of the quarkonia states will be employed for the study of their decay 
rates. 



2.2 Coulomb plus power potential in variational approach 



The heavy-heavy bound state systems such as cc, bb, may also be described by a nonrelativistic Hamiltonian 
given by 0(91 [TO] 

H = M + ^- + V(r), (8) 
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Table 2: Masses of bottomonia in GeV/c 2 



State 


ERHM 




CPP, 




32 


m 


42 


[35] 


[33 








v=0.5 


,=1.0 


v— 1.5 


,=2.0 


PDG 












1^0 


9.425 


9.426 


9.411 


9.399 


9.389 


9.300 


9.457 


9.421 


9.400 


9.454 


9.300 


2 1 S 


10.012 


9.696 


9.826 


9.924 


9.995 




10.018 


10.004 


9.993 




9.974 


3^0 


10.319 


9.824 


10.088 


10.334 


10.529 




10.380 


10.350 


10.328 




10.333 


4^0 


10.572 












10.721 


10.632 








5% 


10.752 












11.059 












9.461 


9.463 


9.468 


9.472 


9.475 


9.460 


9.460 


9.460 


9.460 


9.505 


9.460 


2 3 S! 


10.027 


9.702 


9.841 


9.951 


10.032 


10.023 


10.023 


10.024 


10.023 


10.013 


10.023 


3 3 S 1 


10.329 


9.827 


10.097 


10.334 


10.529 


10.355 


10.385 


10.366 


10.355 


10.335 


10.381 


4 3 S 1 


10.574 










10.579 


10.727 


10.643 




10.577 


10.787 


5 3 S 1 


10.753 










10.865 


11.065 






10.770 


11.278 


1 3 P 


9.839 


9.664 


9.755 


9.820 


9.866 


9.859 


9.894 


9.860 


9.863 


9.855 


9.865 


1 3 P! 


9.873 


9.670 


9.775 


9.820 


9.866 


9.893 


9.941 


9.892 


9.892 


9.875 


9.895 


1 3 P 2 


9.941 


9.683 


9.792 


9.866 


9.913 


9.912 


9.983 


9.910 


9.913 


9.887 


9.919 


l l Pl 


9.907 


9.672 


9.775 


9.852 


9.911 




9.955 


9.900 


9.901 




9.894 


2 3 P 


10.197 


9.803 


10.035 


10.228 


10.379 


10.232 


10.234 


10.231 


10.234 


10.212 


10.238 


2 3 Pi 


10.207 


9.806 


10.044 


10.246 


10.406 


10.255 


10.283 


10.258 


10.255 


10.227 


10.264 


2 3 P 2 


10.227 


9.811 


10.062 


10.282 


10.462 


10.268 


10.326 


10.271 


10.268 


10.237 


10.283 


2 : Pi 


10.217 


9.808 


10.053 


10.264 


10.434 




10.296 


10.263 


10.261 




10.260 
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Table 3: Theoretical predictions of the ground state masses (in GeV) and |i?(0)| 2 of heavy quarkonia 



Models 








M P 


My 


\Rcw(0)\' A 

GeV 3 




Theory 


Expt [31, 


Theory 


Expt [31, 


ERHM [16] 




0.356 


2.985 


2.9804 


3.096 


3.097 


0.556 


BT g] 




0.360 


2.980 


± 


3.097 




0.810 


cc PL(Martin) 


A 


0.430 


2.980 


0.0012 


3.097 




0.999 


Log [5] 




0.370 


2.980 




3.097 




0.815 


Cornell [7] 




0.310 


2.980 




3.097 




1.454 


CPP V 01/ = 


=0.5 


0.300 


3.000 




3.092 




0.610 




1.0 


0.300 


2.950 




3.112 




1.100 




=1.5 


0.300 


2.912 




3.129 




1.508 




:2.0 


0.300 


2.882 




3.144 




1.850 


ERHM [16] 




0.241 


9.452 




9.464 




4.990 


BT g] 




0.241 


9.377 


9.300 


9.464 


9.4603 


6.477 


bb PL(Martin) 


,5 


0.270 


9.398 


± 


9.462 


± 


4.591 


LOG [B] 




0.245 


9.395 


0.002 


9.460 


0.00026 


4.916 


Cornell [7] 




0.217 


9.335 


± 


9.476 




14.05 


CPP V 01/ = 


=0.5 


0.233 


9.426 


0.002 


9.463 




3.908 




1.0 


0.233 


9.411 




9.468 




5.988 




=1.5 


0.233 


9.399 




9.472 




7.728 




=2.0 


0.233 


9.389 




9.475 




9.181 
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where M = tuq + ttiq, M\ — m< J +m q , tuq and tuq are the mass parameters, p is the relative momentum of 
each quark and V(r) is the quark antiquark potential given by [8] 

V(r) = ^+Ar»; a c = \ a s (9) 
r 3 

Here, we use the Coulomb plus power potential where the power index v varies from 0.5 to 2. Each value of v 
corresponds to different potential scheme. Under the variation approach one requires a trial wave function to 
compute the expectation value of the Hamiltonian. In our previous work [8l [9] , we have used the harmonic 
oscillator wave function for the study of the ground state masses, decay constants and decay rates. The 
value of the wave function at origin is very crucial parameter for estimation of decay constant and decay 
rates. But it is found that the radial wave function at r=0 are underestimated with the gaussian trial wave 
function. We also came across problem of overestimation with the Gaussian like wave function (HO basis) 
in the predictions of the orbital excited states. This suggests us to look for an alternate trial wave function 
for the present study. 

Accordingly, we employ hydrogenic trial wave function given by 

*"<" " {^Tw) 1 ' 2 <" r >' '-"" 2 L ™M (10 » 

for the variational calculations in the present study. Here, \i is the variational parameter and L^^^^fxr) 
is Laguerre polynomial. The Potential parameter, A for each choices of v is fixed to yield the ground state 
mass of the meson. 



The expectation value of the Hamiltonian for the ground state is obtained as 

Efa u) = M + - — + - (- mc + A r( " + 3) ) (11) 

The wave function parameter \x is determined here by using virial therom for a chosen value of v. Thus, we 
obtain the spin average mass {Ms a) °f the system from Eqn (jlip . In our earlier study of quarkonium we 
have employed the gaussian trail wave function for the mass spectra, decay constant and decay rates, but 
the predictions of the decay constant and decay rates were underestimated significantly [H [9] . For the same 
region in the present study we are using hydogenic like wave function. Unlike the earlier case of ERHM, 
here we consider the one gluon exchange for the spin-spin and spin-orbit interactions given by [11] 

Vs ^ {r) = l^k-S Q -^ 4n6{r) (12) 

y mQiriQ 

and 

V L -s{r) = q 5- (13) 

6 uiqitiq r a 
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The value of the radial wave function -R(O) for h and 1" states would be different due to their spin 
dependent hypcrfinc interaction. The spin hyperfinc interactions of the heavy flavour mesons are small and 
this can cause very little shift in the value of the wave function at the origin. Many models do not consider 
this contribution to the value of R(0). However, we account this correction to the value of R(0) by considering 



Where (SF)j and (esd)tiJ is the spin factor and spin interaction energy of the meson in the nJ state, 
while R(0) and Ms a correspond to the radial wave function at the zero separation and spin average mass 
respectively of the QQ system. The mass parameters m& = 4.66 GeV, m c =1.31 GeV and the mass difference 
between the pseudoscalar and vector meson due to the chromomagnetic hyperfine interaction is computed as 
described in Ref [§]. The potential parameter A of CPP vA is fixed at 0.276, 0.167, 0.101, 0.060 GeV +1 for 
cc and 0.195, 0.158, 0.126, 0.098 GeV 1/+1 for bb, for the choices of v = 0.5, 1.0, 1.5, 2.0 respectively to get 
the ground state spin average masses of the quarkonia. The mass spectra computed in this model are listed 
in tables [1] and [3 along with other model predictions. The computed ground state spectroscopic parameters 
like the vector and pseudoscalar masses, the square of the radial wave functions at the origin for cc and bb 
systems are tabulated in Table[3]for v = 0.5, 1.0, 1.5 and 2.0. Similar parameters of other potential models 
are also listed. These data will be used to compute the decay rates. 

3 Decay rates of cc and bb mesons in NRQCD formalism 

The decay rates of the heavy quarkonium states, light hadrons, photons or pairs of leptons are among the 
earliest applications of perturbative QCD [55] HS1 HZ] . The decay rates of the mesons are factorized into 
a short-distance part that is related to the annihilation rate of the heavy quark and antiquark and a long- 
distance part containing all nonperturbative effects of the QCD. The short-distance factor calculated in 
terms of the running coupling constant a(M) of QCD is evaluated at the scale of the heavy-quark mass 
M, while the long-distance factor is expressed in terms of the meson's nonrelativistic wave function, or its 
derivatives, evaluated at origin. The di-gamma decay of 1 Sq state and the leptonic decay of 1 state using 
the conventional Van-Royen Weisskopf formula [28] 




Rnj(0) = R(0) l + (SF)j 



{£SD)nJ 
M SA 



(14) 




r 




Rp(0) 



(15) 



and 



Tvw 




R 2 y(0) 



(16) 
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Table 4: + ^ 7 7 and 1 — > l + I decay rates (in keV) of cc and bb mesons. 



AT .-./-I do 




^ — > 7 7 




1 — — 

1 — >■ 






F 

1 o 


1 NRQCD 1 _Ba;pt. ^oJ-J 


1 vw 


1 NB.QCD 


1 Ba;pt. .oJ-J 


H/xvrllVl [Id] 


/ .0( 


o.4Ul 


K ACZCi 

o.4t>y 


O Q77 




BT [4] 


11 1 n 

11.19 


5. loo 7.0_ 9 


8.31 


O T"1 1 

3.711 


r /in 1 nir 

5.40 ± 0.15 


cc (Martin J [5] 


1 o o 1 

13. ol 


a AC A 

O.U54 


y.yo 


4.byi 


±0.0/ 


Log [o] 


11. ZD 


5.004 




O A/IO 

o.y4z 




Cornell [7] 


20.09 


8.614 


14.50 


6.704 




CPP, [5]i/ =0.5 


8.173 


2.213 


6.130 


1.766 




=1.0 


14.649 


5.513 


11.053 


3.017 




=1.5 


19.971 


7.889 


15.165 


3.938 




=2.0 


24.297 


9.920 


18.549 


4.591 




ERHM [12] 


0.440 


0.207 


1.320 


0.620 




BT [4] 


0.569 


0.278 


1.717 


0.833 




66 PL(Martin) [5] 


0.408 


0.195 0.364 [T3] 


1.216 


0.588 


1.314 ± 0.029 


Log [5] 


0.437 


0.211 


1.303 


0.633 




Cornell [7] 


1.258 


0.597 


3.719 


1.801 




CPP„ 8 ^=0.5 


0.345 


0.163 


1.035 


0.489 




=1.0 


0.529 


0.251 


1.587 


0.748 




=1.5 


0.683 


0.324 


2.047 


0.965 




=2.0 


0.811 


0.385 


2.433 


1.146 





Table 5: Leptonic decay widths (in keV) of cc(n 3 Si) 



State 


ERHM 


CPP, y 


32, 


m 


m 


m 


m 


m 






u = 0.5 v=l.O 


PDG 


Pert 


Non-pert 










5.469 


6.130 11.053 


5.55 ± 0.14 


4.28 


1.89 


4.80 


7.82 


6.72 ± 0.49 


V(2 3 ft) 


2.140 


1.232 2.973 


2.48 ± 0.06 


2.25 


1.04 


1.73 


3.83 


2.66 ± 0.19 




0.796 


0.551 1.495 


0.86 ± 0.07 


1.66 


0.77 


0.98 


2.79 


1.45 ± 0.07 


^(4 3 ^i) 


0.288 




0.58 ± 0.07 


1.33 


0.65 


0.51 


2.19 


0.52 ± 0.02 



Table 6: Leptonic decay widths (in keV) of bb(n 3 Si) 



State 


ERHM 


CPP„ 


32 


[39] 


m 




[42] 


02 






v = 0.5 i/ = 1.0 


PDG 








Pert 


Non-pert 


T(l a 5'i) 


1.320 


1.035 1.587 


1.340 ± 0.018 


1.45 ± 0.07 


1.314 




5.30 


1.73 


T(2 3 5*i) 
T(3 3 S 1 ) 


0.628 


0.155 0.390 


0.612 ± 0.011 


0.52 ± 0.02 


0.576 


0.426 


2.95 


1.04 


0.263 


0.066 0.211 


0.443 ± 0.008 


0.35 ± 0.02 


0.476 


0.356 


2.17 


0.81 


T(4 3 5i) 
T(5 3 5i) 


0.104 




0.272 ± 0.029 




0.248 


0.335 


1.67 


0.72 


0.040 




0.310 ± 0.070 




0.310 


0.311 
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as well as using the NRQCD formalism have been computed. The NRQCD factorization expressions for the 
decay rates are given by [13] 



rCSo-77) = 2Imf zf So) l(0|xWSo)| 



2Img~ f ~ f ( 1 So) 



Q 



xRe 



(17) 



TV3o + -\ 2/m/ee( 3 Sl) i, . + ,|3c \ 



2/m.g ee ( 3 S 1 ) 



xRe 
+0(v 4 T) 



8 WW|o)(o|xV(-i£)Vl 3 Si> 



(18) 



The matrix elements that contributes to the decay rates of the S wave states into rj c — » 77 and ^ — > e + e 
through next-to- leading order in w 2 , the vacuum-saturation approximation gives [13j 



( 1 5o|0( 1 S'o)| 1 S'o) = |(0|x t V| 1 S'o)r[l + 0( U 4 r)] 



(19) 



'Siio^soi^i) = |<o|xWl 3 Si>| 2 [i + o(v 4 r)] 



(20) 



( 1 5o|7' 1 ( 1 5 )| 1 5o) = i?e 

+0(v 4 r) 



( 1 5 o |^ x |0)(0|x t (-^) 2 ^r5 o ) 



(21) 



5 5i|P!( 3 5i)| 3 5i) = i? e [( 3 5 1 |^ t ( 7x|0) 

(0| X V(-l5)V| 3 Si> 

+o(w 4 r) 



(22) 



The Vacuum saturation allows the matrix elements of some four fermion operators to be expressed in terms 
of the regularized wave-function parameters given by [13j 



{'SolOCSo^So) = -|ffi So (0)| 2 [l + CV)] 



(23) 
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("SilO^SOfSi) = ^\R 3Sl (0)\ 2 [l + O(v 4 )} (24) 



( 1 5o|^i( 1 5o)r5 ) = -t-\R* cw X7 2 R CW \[1 + 0(v 4 )} (25) 

Z7T 



< 3 ^i|^ 1 ( 3 5' 1 )| 3 5' 1 > = - — \R* CW ^R cw \[l + 0(v 4 )} (26) 

Z7T 

The factorization formula for electromagnetic annihilation, the decay rates for r\ c — > 77 and ^ — > e + e _ are 

r(So^ 77 ) = — —2 — |a So | 



12^»Re(R* V 2 i? s ) 
Trm^ 



+0(w 4 r) (27) 



r( 61 e^e ) = 5 i?3 Sl 

7T 



7V C Im,?^ 3 ^) 



i?e(i?* V 2 i? s ) 



+0(« 4 r) (28) 
The short distance coefficients f 's and g's computed in the order of a 2 as [T3] 

Imf^CSo) = ttQV (29) 
7mp 77 ( 1 S'o) = -^Y~a 2 (30) 
Im/ ee ( 3 S0 = (31) 
Img^Sx) = -^fa 2 (32) 

Where C F = (iV 2 - l)/(2iV c ) = 4/3. 

The decay rates are computed using conventional Van-Royen Weisskopf formula for T & Tyw as well as 
using NRQCD expressions. Here N c = 3 is the numbers of colour and a is the electromagnetic coupling 
constant. The phenomenological mass parameters and the strong running coupling constant employed in 
the meson mass predictions are being employed here to compute the decay rates. For comparison, we have 
listed the results in Tables 0] \5\ and [5] 
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4 Summary and Conclusion 



We have presented the decay rates of the charmonia and bottomonia systems using spectroscopic parameters 
obtained from different potential models within the NRQCD formalism. As the non-perturbative aspect of 
the interaction has already been taken into account in the confinement schemes of quarks and gluons, the 
residual one gluon exchange effects employed in this study are treated perturbatively. The predictions 
of harmonic plus coulomb potential (ERHM)are found to reproduce the experimental results with good 
accuracy. The predicted tj c (2S) mass is in good agreement with recent experimental mass. 

Using the radial wave functions of cc and bb systems of different potential models, the decay rates of ^ — * 77 
and 1 — > are computed within the NRQCD formalism. The results are compared with the values 
obtained using the conventional formula (To, Tyw) &s well as with the respective experimental results. It is 
observed from Table |4] that the computed radiative decay of r/ c from ERHM using the conventional formula 
(To) is closer to the experimental value while all other model predictions using NRQCD formalism found 
to be closer to the experimental values within its error bars on either side. Similar is the case for the pre- 
dictions of the r ee width of J/if>. The leptonic decay widths of higher excited states also found to be in 
good agreement with experimental values. Further, the radiative decay rates of rjb meson predicted by all 
the models except Cornell potential are close to each other. In case of bb, all the models seem to be in good 
agreement with experimental results. 

Though the predictions using conventional formula are far from the experimental results in most cases, the 
predictions based on NRQCD are found to be in accordance with the experimental results. The present 
study in the determination of the S wave masses and decay rates of cc and bb systems provide future scopes 
to study light hadron decays and various transition between the excited states of these mesonic systems 30i . 
It can be concluded that the NRQCD formalism has most of the corrective contributions required for most 
of the potential models for the right predictions of the decay rates. It must also be noted that the ERHM 
model predictions of the decay rates are in excellent agreement with the experimental values. It may be due 
to the fact that its parametric description of the quarkonia through the spectroscopy is more close to reality. 

In case of CPP„, no single value for v is found to account for all decay widths and masses correctly. The 
purpose of varying the potential index v is to observe the variation of the different properties of the quarkonia 
with the choice of where different potential forms. From the present study it is also found that the potential 
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index v is not the only factor affecting estimation of masses and decay rates as the wave function at r = 
is an important parameter. The values of the decay rates are overestimated with hydrogenic wave function 
while largely underestimated with gaussian trial wave function [9]. Hence, the present study indicates that 
the right form of the wave function lies between the hydrogenic to gaussian i.e. of the form e~ MrP with 
1 <p < 2. 

The present study also suggests that the spectroscopic parameters perhaps are not sufficient for understand- 
ing the decay dynamics of hadrons. 
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